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Abstract 

The physical relevance of the fractional time derivative in quantum mechanics 
is discussed. It is shown that the introduction of the fractional time Scrodinger 
equation (FTSE) in quantum mechanics by analogy with the fractional diffusion 
H ~^ '^^^ lead to an essential deficiency in the quantum mechanical descrip- 
tion, and needs special care. To shed light on this situation, a quantum comb 
model is introduced. It is shown that for a = 1/2, the FTSE is a particular 
case of the quantum comb model. This exact example shows that the FTSE is 
insufficient to describe a quantum process, and the appearance of the fractional 
time derivative by a simple change ^ — )■ in the Schrodinger equation leads 
to the loss of most of the information about quantum dynamics. 
PACS: 05.40.-a, 05.45.Mt 
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1. Introduction 

In quantum physics, the fractional concept can be introduced by means of 
the Feynnian propagator for non-relativistic quantum mechanics as for Brown- 
ian path integrals 1,]. Equivalence between the Wiener and the Feynman path 
integrals, established by Kac Q], indicates some relation between the classical 
diffusion equation and the Schrodinger equation. Therefore, the appearance of 
the space fractional derivatives in the Schrodinger equation is natural, since both 
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the standard Schrodinger equation and the space fractional one obey the Markov 
process. As shown in the seminal papers |3|, Ij], it relates to the path integrals 
approach. As a result of this, the path integral approach for Levy stable pro- 
cesses, leading to the fractional diffusion equation, can be extended to a quan- 
tum Feynman-Levy measure which leads to the space fractional Schrodinger 
equation 0,0]. 

The fractional time Schrodinger equation (FTSE) was first considered in Q , 
where a fractional time derivative was introduced in the quantum mechanics by 
analogy with the fractional Fokker-Planck equation (FFFE), by means of the 
Wick rotation of time t — > —it/h . The dynamics does not correspond to the 
unitary transformation: the Green function is found in the form of the Mittag- 
LefHer function and does not satisfy Stone's theorem on one-parameter unitary 
groups A generalization for the space-time fractional quantum dynamics 
[3, 0| was performed and a relation to the fractional uncertainty [9] was studied 
as well. It was also shown that the FTSE introduces new nonlinear phenomena 
in the semiclassical limit, and this semiclassical approach differs from those 
described in the framework of the standard Schrodinger equation [l^ . 

The fractional time quantum dynamics with the Hamiltonian H{x) is de- 
scribed by the FTSE 

{fhr^^^^-m^,t), (1) 

where a < 1. For agreement of the dimension in Eq. (HI), all variables and 
parameters are considered dimensionless, and h is the dimensionless Planck 
constant, see also . For a = 1, Eq. ([T]) is the "conventional" (standard) 
Schrodinger equation. For a < 1 the fractional derivative is a formal notation 
of an integral with a power law memory kernel of the form 

d'-m ^ ji-c dm _ r (^-t)-" om 

dt- - * dt ~ Jo r(l-a) dr ' 

which is the Caputo fractional derivative: ^ g^i*-* = QD^ip{t) (see Appendix). 

In this paper we present an exact example where the fractional time deriva- 
tive is naturally introduced and has a well defined physical meaning. We con- 
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sider a case with a = 1/2, when fractional quantum dynamics can be modelled 
by means of the conventional quantum mechanics in the framework of a comb 
model, and the FFPE for a = 1/2 is a particular case of the quantum comb 
model which is a quantum counterpart of a diffusive comb model . We 

also show that the FTSE ([2]) is insufficient to describe a quantum process, and 
the appearance of the fractional time derivative leads to the loss of most of the 
information about quantum dynamics. The main idea is to show that Eq. ([1]) 
is a result of a projection of the two-dimensional (x, y) comb dynamics in the 
one-dimensional configuration space. For diffusion, which is described by the 
Fokker-Planck equation, this projection is a simple integration over the y space. 
In quantum mechanics this projection is performed by means of the Fourier 
transform 'i>{x,y) — > ^i{x) = J-y'i>{x,y), where / is the Fourier index. For 
the comb model this procedure can be treated exactly, and we shall show that 
Ec^. ([l} is valid only for the zero Fourier component 'ijj{x) = ^q{x). All other 
components are not described by the FTSE ([1]) and lost in the framework of 
this equation. The diffusive comb model is an analogue of a Id medium where 
fractional diffusion has been observed [ll[ [3] ■ It is a particular example of a 
non-Markovian phenomenon, explained in the framework of the so-called contin- 
uous time random walks (CTRW) QQQ. This model is also known as a toy 
model for a porous medium used for exploring of low dimensional percolation 
clusters |15|. 



2. Quantum Comb Model and FTSE 

A special quantum behavior of a particle on the comb can be defined as the 
quantum motion in the d+1 configuration space (x, y), such that the dynamics in 
the d dimensional configuration space x is possible only at y = 0, and motions in 
the x and y directions commute. Therefore the quantum dynamics is described 
by the following Schrodinger equation 

z/.-=%Mx)vI.--_, (3) 



3 



where the Hamiltonian H = H{'x) = — + V(x) can be different from clhH 
in Eq. ([l}. It governs the dynamics with a potential V^(x) in the x space, while 
the y coordinate corresponds to the Id free motion. All the parameters and 
variables are dimensionless 0- 

By analogy with the diffusion (classical) comb model, we are concerned with 
the dynamics in the x space. But simple integration of the wave function over 
the y coordinate is not valid. Therefore, one carries out the Fourier transform 
in the y space J\,5'(x, = '^(x, = and as a result of this, Eq. ([3]) 
reads 

z/i^=i?(x)*(x,0,t) + ^*,. (4) 

To obtain this equation in a closed form, one needs to express the wave function 
at y = \I'(x, 0,t) by the Fourier image To this end the Laplace transform 
of Eq. is performed £[5'(x, — ^s{^,y)- The solution in the Laplace 
domain reads 

*,(x,j/) = *,(x,0)exp + , (5) 

where we used ^/2i — {1 + i). Performing the Fourier transform ^s.i{^) = 
T['^s{^,y)] one obtains from Eq. ([5]) 

, ^ rrr, , 2i(l + i)\/s/Ji _ 
vl/,,(x) = vI;,(x,0)J-,e^(i+^)V«/M.I = V ^^(^^p) (g) 

— 2is/h 



Finally, the Laplace inversion for 5's(x, 0) determines the wave function at j/ — 

^ P-2is/h - 
2i{l + i)\ls/}i 



vl/(x,0,i) = ^.,(x) '^'''tA ■ (7) 



1 Analogously to the FTSE {TJ, following Ref. fl, one introduces the Planck length Lp 



y'hG/c^, time Tp = ^ YiG / c' , mass Mp = ^jTiclG^ and energy Ep = Mpc^, where ti, G, 
and c are the Planck constant, the gravitational constant and the speed of light, respectively. 
Therefore, quantum mechanics of a particle with mass m is described by the dimensionless 
units x/Lp X, y/Lp — > y, t/Tp — > t, while the dimensionless Planck constant is defined 
as the inverse dimensionless mass h = Mp/m. Note, that the dimensionless potential is now 
y{x) ^ V{x)/Mpc2. 
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Let us, first, consider a simple case with I = 0. We have from Eq. ([7]) 
^s(x, 0) — ^ \J s,a- Then we define ^'o(x,f) — ip^x), and, carrying out 
the Laplace transform in Eq. (j?]), we arrive at the definition of the Caputo 

1/2 

fractional derivative ^ in the Laplace domain (see Appendix) C[aDf' ^{t)] = 
s^/^^/j{s) - s-i/2V'(0). Finally, carrying out the inverse Laplace transform and 
redefining — ^ H , one obtains the FTSE which coincides exactly with Eq. 
O for a = 1/2. 

Repeating this procedure for an arbitrary I, one performs the Laplace trans- 
form of the term proportional to l'^ j %fs in Eq. Performing simple operations 
of fractional calculus and taking into account Eqs. (|A. 3p and (jA. 9[) . one ob- 
tains 

This comb FTSE describes the quantum dynamics in the x configuration space. 
The index I corresponds to an effective interaction of a quantum system with an 
additional degree of freedom, while the fractional time derivatives, with a = 1/2, 
reflect this interaction in the form of non-Markov memory effects. 



3. Green's Function 

The initial value problem with the initial condition ^(t = 0) = ^'o(x, ?/) 
is described by Green's function. For the complete analogy with the classical 
comb model and fractional diffusion [l^ , the boundary conditions for the y 
direction defined at infinities y = ±oo are ^'(x, i) = 9^11 (x, t)/dy = 0. Using the 
eigenvalue problem 

i?(x)7AA(x)-AV'A(x), (9) 

we present the wave function in Eq. ^ as the expansion ^'(x, y, t) — (t^xiy, ^)V-'a(x), 
where also supposes integration on A for the continuous spectrum. For the 
fixed A we arrive at the dynamics of a particle in the S potential 
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Taking into account that the Green function of Eq. ([3]) has the spectral decom- 
position, 

G(x, y, t- x'y') - ^^^2^' ^5 2/')V'I(x)V^a(x') , 

A 

we obtain that the Schrodinger equation for the Green function G\{y, t; y') reads 
- dGx d^Gx 



ih- 



dt 



2 dy^ 



X6iy)Gx + ih5{y)5{t) . 



(11) 



Here the initial condition is already taken into account. The Green function 
for this Schrodinger equation has been obtained in [l^, [3| , for free boundary 
conditions at infinities. For the chosen boundary conditions it is instructive to 
employ Eq. ([5]) in the Laplace domain. Then replacing the eigenvalues A by the 
Hamiltonian -ff(x), one obtains the Green's function in the form of the inverse 
Laplace transform 



G 



H{x),y,t) 



£-1 



(12) 



-H{yi)~i{l + i)Vh3s' 
One performs the inverse Laplace transform, using the following presentation 
for the denominator 



exp| -u[i?(x) -i(l + i)yPI]| 



du . 



This presentation is valid for any spectrum A due to the second term in the 
exponential. Therefore, the Green function reads 



G 



ff(x),y,i) 



ih 



{\y\ + h^u)exp -m7J(x)- 



i{\y\ + h^uf 



V2^ 

Using the Fourier transform for the exponential 



2ht 



du. (13) 



exp 



. ii\y\ + h^ur 
2ht 




^thte/2-t(,{\v\+h''u) 



one presents the Green function in the following convenient form 



G(x,y,t:x') — - — 
^ ' ^ 27rt 



exp 



Mff (x) - ie,{\y\ + h^u) du}d^. (14) 



This expression is convenient for further analysis in the framework of the path 
integral. 
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4. Path Integral Presentation 



As an example, it is instructive to consider a free particle, because it has a 
straightforward relation to the original diffusive comb model [lllll^l- Therefore, 
we find the Green function along the structure x axis in the coordinate space 
G{x,y,t;x') — {x'\G H{x.),y,t \x) for a free particle of a unit mass with the 
Hamiltonian H = p'^/2. Expressing the exponential {x'\e~'"^\x) in the path 



integral form, 

G{x,y,t;x') = 



, "'"^ e 2uh^ J one obtains from Eq. (ITil) : 

V 2Trh^u 



01) r 



exp 



ix-x'Y' 
2h?u 



Integration over complex "time" u yields the following expression 



I{A, B) 



exp[ ouj— ;= 



(15) 



(16) 



where A = ^ ^^'^ ^ — -~iS,h^- Differentiation of Eq. ((TE)) with respect to 
A and B yields the equation ^ oBdA^ ^ '^)- ^ solution of this equation is 

I{A,B) 



: exp[2V AB] 



4VAB 

Performing this integration, we arrive at integration over ^ that is carried out 
in the stationary phase approximation for the long time asymptotics ht ^ 1. 
This yields 

4MJ 



^^(Co)exp 



(17) 



where the stationary point is ^0 — Taking integration (jl7p into account, 
we finally, obtain the Green function in the form 

,•1/4 



G{x,y,t;x') 



4:Trty^2\x - x'\ 
X exp 



_,,d)-%,„,-,-.W|i 



i— — h i\/ i\v\/ht\x — x'\ 
2ht ^ ' 



(18) 



This solution also satisfies the boundary conditions at a;, y = ±00, where the 
Green function vanishes. 
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5. Conclusion 



Physical relevance of the fractional time derivative in quantum mechanics is 
discussed. It is shown that the introduction of the fractional time Scrodinger 
equation in quantum mechanics by the change ^ — >■ by an analogy with the 
fractional diffusion can lead to essential deficiency of the quantum mechanical 
description, and needs special care. To shed light on this situation, a quantum 
comb model is introduced. We observed that the fractional time derivative, at 
least for a — 1/2, reflects an effective interaction of a quantum system with an 
additional degree of freedom. In the classical diffusion comb model, diffusion 
in the y direction is responsible for traps that lead to subdiffusion along the x 
structure axis, and this phenomenon is described by the time fractional deriva- 
tive This description in the framework of the fractional Fokker-Planck 

equation is identical to the diffusion comb model |12l . I16| . In the quantum case 
the situation differs essentially from fractional diffusion. First of all, the quan- 
tum comb model ([3|) and the FTSE ([T]) are not identical. As shown here, the 
FTSE is an equation only for the zero Fourier component of the wave function, 
and this equation is insufficient to describe the complete dynamics in the x space 
of the system. The FTSE is a particular and limiting case of the comb fractional 
time Schrodinger equation ([8]), obtained here, which describes the quantum dy- 
namics in the x configuration space, and the Fourier index I corresponds to 
an effective interaction of the quantum Hamiltonian -ff(x) with an additional 
degree of freedom, while the fractional time derivatives, with a — 1/2, reflect 
this interaction in the form of non-Markov memory effects. This equation can 
be readily solved by the Laplace transform and using the eigenvalue problem 
H{x)^\{x) = Xip\{x.). Nevertheless, we obtained the Green function directly 
from the quantum comb model of Eq. ([3]) in a form suitable for the paths in- 
tegral presentation. An example of H{'x) that corresponds to a free particle is 
considered. Note also that this expression for the Green function is also suitable 
for the semiclassical treatment of more complicated Hamiltonian systems -ff(x). 

In conclusion, we admit that this exact example shows that the FTSE ^ is 
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insufficient to describe a quantum process, and the appearance of the fractional 
time derivative by a simple change ^ — in the Schrodinger equation leads 
to loss of most of the information about quantum dynamics. 

This work was supported in part by the Israel Science Foundation (ISF) and 
by the US-Israel Binational Science Foundation (BSF). 



Appendix: Fractional Calculus Tools 

Fractional derivation was developed as a generalization of integer order 
derivatives and is defined as the inverse operation to the fractional integral. 
Fractional integration of the order of a is defined by the operator (see e.g., 



21|) 

aKJix) - -i- / j{v){x - yy-'dy , (A. 1) 



r(a) Ja 

where a > 0, x > a and V(z) is the Gamma function. Therefore, the fractional 
derivative is defined as the inverse operator to a/", namely aDxfi^) — alx'^fi^) 
and alx = aD~". Its explicit form is 

aD^fix) = r fivKx ^ y)-'-^dy . (A. 2) 



r(-a) Ja 

For arbitrary a > this integral diverges, and as a result of this a regularization 
procedure is introduced with two alternative definitions of aD'^. For an integer 
n defined as n — 1 < a < n, one obtains the Riemann-Liouville fractional 
derivative of the form 

jn 

aDlJ{x) = —J-~-f{x) , (A. 3) 



and fractional derivative in the Caputo form (see also 22| ) 



aD^^fix) = airV^"^ {X) , /(") (X) ^ ^/(^) • (A. 4) 

There is no constraint on the lower limit a. For example, when a = 0, one has 

^^-°r(/3 + i) 

ol^RL^ - r(/3 + l-a) ■ 
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We also have 



x'^~°' 



and ar>g[l] = 0, while qD'^lI'^] = x-"/r{l - a). When a = -oo, the resulting 
Weyl derivative is 

-ooVV" = -aoD^ = -oo-D^i = -oqDq . (A. 5) 

One also has _oo-D^e^ = e^. This property is convenient for the Fourier trans- 
form 

= (iA;r/(fc), (A. 6) 

where T[f{x)] = f{k). This fractional derivation with the fixed low limit is 
also called the left fractional derivative. However, one can introduce the right 
fractional derivative, where the upper limit a is fixed and a > x. For example, 
the right Weyl derivative is 

The Laplace transform of the Caputo fractional derivative yields 

CioD^fix)] = s"f{s) ~J2f^'^ (0+)s«-i-'= , (A. 8) 

fe=0 

where C[f{x)] = f{s), which is convenient for the present analysis, where the 
initial conditions are imposed in terms of integer derivatives. We also use here 
a convolution rule for < a < 1 

C[I^f{x)]=s-y~{s). (A. 9) 
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